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Abstract 

Hydrogen reduction of lunar regolith has been proposed as a viable technology for oxygen production 
on the moon. Hydrogen reduces FeO present in the lunar regolith to form metallic iron and water. The 
water may be electrolyzed to recycle the hydrogen and produce oxygen. Depending upon the regolith 
composition, FeO may be bound to Ti02 as ilmenite or it may be dispersed in glassy substrates. Some 
testing of hydrogen reduction has been conducted with Apollo-returned lunar regolith samples. However, 
due to the restricted amount of lunar material available for testing, detailed understanding and modeling 
of the reduction process in regolith have not yet been developed. As a step in this direction, hydrogen 
reduction studies have been carried out in more detail with lunar regolith simulants such as JSC-1 A by 
NASA and other organizations. While JSC-1 A has some similarities with lunar regolith, it does not 
duplicate the wide variety of regolith types on the moon, for example, it contains almost no ilmenite. 
Nonetheless, it is a good starting point for developing an understanding of the hydrogen reduction process 
with regolith-like material. In this paper, a model utilizing a shrinking core formulation coupled with the 
reactor flow is described and validated against experimental data on hydrogen reduction of JSC-1 A. 

Nomenclature 

c molar concentration of hydrogen in particle (moles/m 3 ) 

Ci molar concentration of water (moles/m 3 ) 

c s molar concentration of hydrogen at reaction surface (moles/m 3 ) 

c x molar concentration of hydrogen at location x in reactor (moles/m 3 ) 

c 0 molar concentration of hydrogen at upstream end of regolith bed (moles/m 3 ) 

D effective gas diffusion coefficient (m 2 /s) 

AH endothermic heat of reaction (Joules/mole) 

k equilibrium constant 

L length of regolith bed (m) 

N number of regolith particles in the bed 

P gas pressure (Pa) 

r regolith particle radial coordinate (m) 

rp regolith particle radius (m) 

r s radius of reaction surface (m) 

R r internal radius of reactor (m) 

t time (s) 

tj) intra-particle diffusion time (s) 
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tR 

reactor residence time (s) 

T 

temperature (K) 

u 

velocity through reactor (m/s) 

F bed 

volume of regolith bed (m 3 ) 

Vvoid 

void volume in regolith bed (m 3 ) 

X 

axial coordinate (m) 

a 

particle conversion 

e 

void fraction in regolith bed 

Pa 

molar concentration of FeO in regolith (mol/m 3 ) 


Introduction 

A major objective of lunar in-situ resource utilization (ISRU) (Refs. 1 to 3) is production of oxygen 
from lunar regolith. Mare basalts containing oxygen include olivine, ilmenite, and pyroxene (Ref. 4). 
Energy has to be supplied to remove the oxygen, i.e., the associated chemical reactions are endothermic. 
Free energy considerations show that the iron containing minerals require the least energy input for this 
purpose (Ref. 5). Various technologies have been proposed for the chemical conversion. These include 
ilmenite reduction by hydrogen, molten electrolysis, and carbothermal reduction using methane (Refs. 6 
and 7). The reduction processes involve solid-gas reactions. Such systems are well developed for earth 
applications but mainly for the purpose of metal extraction, as opposed to water vapor and oxygen 
extraction for lunar applications. In addition, in the lunar context, the reactor must interface with the other 
sub-system processes such as upstream regolith extraction and beneficiation and downstream electrolysis 
and phase separation in a way that establishes the most favorable balance between efficiency, robustness, 
and equivalent system mass. A need exists, therefore, to establish an analytical framework to couple these 
sub-systems together. The building blocks of this framework are the individual sub-system models. This 
paper discusses the development and validation of a generalized analytical model of the chemical reactor 
sub-system. 

While the developed methodology should be applicable to different gas-solid reactors, only the 
reduction of FeO by hydrogen will be considered here. FeO occurs in lunar regolith as ilmenite, 
particularly in titanium-rich mare basalts, or may be found dispersed in the glassy phase of lunar regolith 
particles. The overall gas-solid reaction considered between FeO and hydrogen is given by: 

FeO + F ^ Fe + 1 1 jO 

The formed water is electrolyzed in the electrolyzer subsystem to produce oxygen. 

The temperature of the reaction is one of the important parameters of the process. The sintering and 
melting temperature of lunar regolith is noted to be in the range of 1370 to 1670 K (Refs. 6 and 8). 
However, hydrogen reduction has been found to be effective at lower temperatures also, for example, in 
the range of 970 to 1270 K (Ref. 9). Therefore, while three-phase (gas-liquid-solid) reactors for hydrogen 
reduction may be applicable at post-melt temperatures, these are not considered here and attention is 
restricted to gas-solid interactions. 

The equilibrium characteristics of the above reduction reaction have been studied to some extent. An 
expression for the equilibrium constant as a function of the reaction temperature alone has been provided 
(Ref. 10): 


k = i 0 (- 21261 ^+°- 6439 ) 


( 1 ) 
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More recently, computational studies have suggested a pressure dependence of the equilibrium 
constant (Ref. 11) but the pressure sensitivity has been found to be much less than the temperature 
sensitivity. This is reasonable since gas- phase moles are conserved in the reaction. The reaction rate and 
chemical kinetics of the reduction process appear to be much less studied. There is consensus that the 
reaction rate is of Arrhenius type and an activation energy for the forward reaction has been estimated as 
18.8 kcal/mole (Ref. 12). As noted earlier, the conversion is endothermic and energy needed for 
conversion, A H, has been estimated to be 1 1 kJ/mol (Ref. 9). 

According to Equation (1), the equilibrium constant for the temperature range of interest is of the 
order of 0.1. This suggests that substantial excess hydrogen beyond the stoichiometric value is needed to 
completely reduce the regolith. Therefore, a closed system where a batch of regolith and a fixed amount 
of hydrogen are placed in a vessel and allowed to react is not practical because of the large volume of 
hydrogen required to provide the needed oxygen for realistic regolith batch amounts. Thus, an open 
system where a batch of regolith is treated with flowing hydrogen is considered here. 

Some hydrogen reduction testing of Apollo-returned samples have been conducted. However, the 
amount of lunar regolith available for testing and analysis is restricted so that a detailed understanding of 
the various processes involved has not been developed. To get around this restriction, NASA has 
developed different regolith simulants from earth-based soils. A comparison of available simulants with 
some of the Apollo-returned samples has been provided (Ref. 13). Recently, NASA-sponsored hydrogen 
reduction studies have been conducted with the simulant JSC-1 A. JSC-1 A has an FeO content of 
1 1.2 percent by weight of simulant which would result in 2.8 percent water upon complete conversion 
with hydrogen (Ref. 13). 

The primary objective of this paper is to develop and validate a model of the hydrogen reduction 
process. A model was previously formulated by the authors (Ref. 14) in terms of the following overall 
steps: (1) inflow of pure hydrogen into the regolith containing reactor, (2) reaction of hydrogen with the 
FeO portion of the regolith, and (3) outflow of left-over hydrogen and produced water vapor. Step (2) is 
the gas-solid reaction step which is the key to the conversion. Results of regolith particulate analysis 
conducted on Apollo-returned samples (Ref. 6) suggest that the particles are not completely solid but have 
some porosity that would enable gases to diffuse into the particles, i.e., the FeO portion can be essentially 
completely converted with the reduced iron remaining behind in the particulates. It may also be noted that 
JSC-1 A has a large glass content with FeO dispersed in the glass phase. The diffusion of gases in the 
glass phase is expected to be the controlling feature for the reduction process. These observations support 
utilization of a “shrinking core” diffusion-reaction model (Ref. 15) for the gas-solid reduction process. 
This model is then coupled with the inflow/outflow of gases to complete the reactor sub-system model. 

In this paper, the model is considered first and a closed-form expression for the FeO conversion is 
obtained for the case of a well-stirred reactor. Two key parameters appear in the expression: (1) a 
diffusion time related to the diffusion of gases in the interior of the regolith particles, and (2) a residence 
time of the gases in the reactor. It is found that for high levels of conversion or for a relatively small ratio 
of residence time to diffusion time, the governing parameter is the diffusion time. Validation experiments 
have been carried out at Pacific Northwest National Laboratory (PNNL) (Ref. 16) and at Lockheed 
Martin (LM) (Ref. 1 7) and are briefly described in Section 3. A factor F defined in Section 2 which is 
related to the diffusion times for different particle sizes is derived from the PNNL data in conjunction 
with the model. Following this step, the obtained factor F is input to the model and used to predict the 
conversion rates in the LM experiments. Although the LM experiments are at significantly larger scale 
than the PNNL experiments, the comparison of model results with the LM experiments is remarkably good. 

Model Description 

The shrinking core model as envisioned here is schematically depicted in Figure 1. Conversion 
proceeds from the exterior particle surface inward. As the reaction progresses into the interior of the 
particle, the reactant and product gases must diffuse through the particle matrix structure formed after the 
conversion. 
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Unreacted core 


Radius of reaction surface, r s (t) 


Particle radius 



Figure 1. — Shrinking core model configuration. 

The diffusion of hydrogen in the particle interior between r p and r s is considered. Assuming quasi- 
steady behavior, and an effective diffusion coefficient, D , for the gases in the interior particle matrix, the 
concentration, c, of hydrogen varies according to 


— (r 2 D— ) = 0 
dr dr 


( 2 ) 


Boundary conditions for c have to be provided at the core (reaction) surface radius, r s , and at the 
exterior particle radius, r p . Note that the particle typically preserves its original (unreacted) size even 
during the reaction process as the main constituents (i.e., the non-FeO portion) are largely unaffected. 

The hydrogen flux is related to the conversion of FeO at the reaction surface. With p a being the molar 
concentration of FeO in the unreacted particle, the conversion of FeO by the flux of hydrogen is given by 


d 4 3 2 

p a nr =4nrD—] 

dt 3 a s 5 dr s 


(3) 


Solving for r s yields 


3[1 - a ] 2/3 - 2[1 - a] - 1 


6D[c x -c s ]t 

Pa r p 


(4) 


where a is defined as 



(5) 


For an inwardly directed diffusive flux of hydrogen to exist within the reacted portion of the particle 
matrix, the value of [ c x - cj must be greater than zero. Assuming that the rate of gas diffusion in the 
particles is much slower than the reaction rate, hydrogen and water may be taken to be in equilibrium at 
the reaction surface resulting in the following relationship (Ref. 14): 


(1 + *) 


(6) 
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It is assumed that the regolith particles form a bed in the reactor. The packing of the particles is 
characterized by the linear loading, NIL, of the particles. For example, the loading for a packed bed will 
be greater than that for a fluidized bed. It is assumed that the gas flow through the bed is characterized by 
a bulk flow velocity, u and the species (i.e., hydrogen or water) concentration is a function of only the 
distance along the bed length, x. The hydrogen concentration, c x , is a function of x due to depletion by the 
reaction. Neglecting inter-particle gas diffusion in the bulk flow, the convection of hydrogen is related to 
its depletion by reaction in the following manner: 


n 2 dc x 2 N 

rt r u- = 4*r,D(-) r - 


(7) 


In the well-stirred reactor, both the gas species concentrations and a are independent of location. The 
value of c x may then be identified with c(L) and the solution of Equation (7) in conjunction with 
Equations (5) and (6) yields: 


c(L) = -\{ 1 + 


k 


l + £ n an 1 -a] m L. 

[1 + 3(1 -£)(—) 1 \i /3l - 

r; [l-o -a) ]« 


( 8 ) 


and, correspondingly, the water concentration at the outlet is given by 

c,(i) = — {1 

' + * P+ 30-«x4) [1 -“ ] 


1/3 


-} 


>r[i-a -cc) m ]u 




(9) 


where 8 is the bed void fraction. The equation for a is obtained by substituting Equation (8) for c x in 
Equation (4): 


3[1 — ot] 2/3 — 2[1 — ot] — 1 = — 


Ft 


[l + 3(l-8)(4) 


[1-oc] 


1/3 


r~ [1 — (1 — ot) 1/3 ] W 


where, the factor F is given by 


F = 6— 

Pa 


k D 
(l + k)~r^ 


( 10 ) 


( 11 ) 


For the case when the residence time, Liu, of the gases is small compared to the diffusion time, r p 2 /D, 
of gases in the particle matrix the equation for a becomes 


3[1 - a] 2/3 — 2[1 — ot] — 1 = -Ft 


( 12 ) 


It may be noted that the time required for high conversion (i.e., as a approaches unity) is also given 
by Equation (12). In practice, this means that unless the residence time is large compared to the diffusion 
time, the time for complete conversion is governed only by the diffusion time of the gases inside the 
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particles. This implies that the factor F is a characteristic descriptor of the reduction process. Examination 
of the expression for F allows for the following observations: 


(a) F has units of s _1 

(b) The equilibrium constant is for Reaction (A) which is, presumably, not simulant dependent. 

(c) The FeO concentration, p OJ is simulant dependent but should be relatively invariant from sample 
to sample of the same stimulant. 

(d) The hydrogen concentration, c 0 , in the reactor is a parameter of the experiment and is operator 
dependent; it is a function of the reactor pressure, P, and temperature, T. 

(e) The particle size, r p , is an effective particle radius. There is generally no ambiguity about particle 
size for a relatively narrow size range; however, when the particle size range is broad, experiments need 
to be done to characterize the effective particle size. One way of doing this is to compare the conversion 
time for a well-characterized particle size with that of the size sample of interest as will be noted in the 
Validation Experiments Section below. 

(f) If the factors in (a) to (d) above are known, then F may be directly related to D, the effective 
diffusion coefficient of the gases in the particle by Equation (1 1). 


Validation Experiments 

Results from the developed model have been compared with data provided by NASA. Two sets of 
experiments were considered: 

(1) The first set of experiments was conducted at PNNL (Ref. 16) in a tube reactor with an internal 
diameter of 1 .64 cm and placed in an oven heated to the desired temperature. The simulant batch size was 
5 grams. Argon was flowed through the sample during the heating process and served to remove any 
trapped moisture in the sample. Upon reaching the desired temperature, the flow of argon was turned off, 
and the hydrogen flow was established. During the reduction process both the hydrogen flow and 
temperature were kept constant. Experiments were conducted for two flow rates of hydrogen (equivalent 
to 0.5 and 1 cm/s velocity under standard conditions), and for reduction temperatures in the range of 800 
to 1000 °C. The reduction time was approximately 4 hr in most cases. In order to investigate JSC-1 A 
particle size effects, tests were conducted with particles sieved for three different sieve sizes as well as 
with unsieved JSC-1 A. The test matrix is shown in Table I where the test points are numbered as in the 
original PNNL report for traceability. 

The PNNL test results were used primarily to assess whether the model captured the reduction 
process and water production rate characteristics and to obtain quantitative values for the factor F. 


TABLE I.— PNNL TEST MATRIX 
[400 mesh = 38 gm, 200 mesh = 75 gm, 115 mesh = 106 gm, 
and 60 mesh = 250 gm] 


Run no. 

Sorting 

Temperature, 

°C 

Gas velocity, 
std cm/s 

74 

None 

1000 

1.0 

76 

None 

900 

1.0 

77 

None 

800 

1.0 

80 

None 

900 

0.5 

81 

None 

1000 

.5 

82 

None 

800 

.5 

87 

>400 mesh 

1000 

0.5 

88 

>400 mesh 

1000 

0.5 

90 

200 to 400 mesh 

1000 

0.5 

91 

60 to 115 mesh 

1000 

0.5 
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Figure 2. — Measured and model predicted conversion versus 
time for JSC-1 A sub-400 mesh particles at 1000 °C. 

(2) The second set of experiments were conducted by LM in a rotating reactor configuration 
(Ref. 17). The reactor section was cylindrical with an inner diameter of approximately 15 cm and length 
of approximately 25 cm. The simulant (unsieved JSC-1 A) loading was 1.5 kg/batch. The reactor and 
simulant were brought up to temperature by external heating. Gas flow was provided through a central 
tube. The inlet flow was impacted onto a plate and suitable placement of the exhaust tube in conjunction 
with reactor rotation appeared to mix the gas well with the regolith. The rotation rate of the reactor was a 
few rotations per second. Argon was flowed during the heat up to remove any trapped moisture. 
Subsequently, the inert gas flow was turned off and the hydrogen flow (3 to 5 slpm) was turned on. 

The LM test results were used to validate the model quantitatively using the value of F obtained from 
the PNNL experiment set. 


Results and Discussion 

PNNL Data Analysis 

The PNNL data in conjunction with the model are utilized to investigate the effects of particle size, 
temperature, and velocity and obtain quantitative values for the factor F (see Eqs. (11) and (12)). 

Baseline Results 

First consider the data from Runs 87 and 88 of the PNNL matrix as a baseline. The two test points are 
essentially the same given the statistical nature of the particle sizes and shapes and serve to verify the 
repeatability of the results. In addition, these tests utilize the smallest particle sizes and the highest 
temperature (1000 °C) in the test matrix and are expected to take the least time for complete conversion. 
In Figure 2, a conversion of unity represents a water production of 2.8 percent by mass of the simulant 
sample, as noted earlier. These small particles are all reduced completely in the reaction time of about 
4 hr. Equation (12) is used to find a fit for F that recovers the behavior of the conversion versus time 
curve during the latter stages of the reduction process. A good fit for the two sets of data is given by 
F= (1/15,000) s~ l . Using this value of F, the model prediction for the conversion is also plotted as a 
function of time and compared to the experimental results for Runs 87 and 88 in Figure 2. 

It is seen that with this value of F, Equation (12) somewhat overpredicts the conversion initially but 
the conversion after a time of approximately 45 min is captured very well. For the small particle sizes 
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used here, the first 45 min represents a conversion of about 0.6. The difference between the model and 
tests 87 and 88 is less than 5 percent at later times. 

Effects of Particle Size 

Tests 88 (<38 pm sieve), 90 (38 to 75 pm sieve), and 91 (106 to 250 pm sieve) provide information 
on conversion times at the same temperature (1000 °C) for different particle sizes (i.e., effective particle 
diameters). The experimentally obtained values of conversion versus time for the different sizes are 
plotted in Figure 3. As is clearly seen, the conversion at a given time is less as the particle size increases. 
The values of the factor F for the larger sizes (Runs 90 and 91) may be obtained from fitting of the 
conversion versus time curve as before (e.g., see Fig. 4 for the case of Run 91) and are found to be: 

F= (1/58,400) s" 1 for Run 90 
(1/175,550) s" 1 for Run 91 



Figure 3. — Measured conversion versus time for JSC-1 A 
particles of different sizes at 1000 °C. 



Figure 4. — Measured and model predicted conversion versus time 
for JSC-1 A 60 to 1 15 mesh particles at 1000 °C. 
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Other than the particle size, the test conditions are identical. Thus, the above variation in F is caused 
by changes in the particle size alone. The effective particle sizes for Run 90 and 91 may be estimated 
from the experimentally derived values of F and from Equation (1 1) by noting that the only variable in 
these tests is the particle radius r p . It is found that for Run 90, the effective particle size is 2.0 times the 
effective particle size for Run 88 (and Run 87) and for Run 91 it is 3.4 times the effective particle size for 
Run 88. For example if we let the effective particle size for Run 88 be 38 pm ( the largest value allowed 
by the corresponding mesh), then, according to the model, for Run 90 it is 76 pm and for Run 91 it is 
129 pm. The value of 76 pm for Run 90 is practically coincident with the upper size for the mesh. 
However, the value of 129 pm is smaller than the 250 pm upper size for the corresponding mesh. It is 
possible that the larger particle sizes in the simulant deviate appreciably from spherical allowing for a 
smaller characteristic size. 

Results for Unsieved JSC-1 A at 1000 °C 

Run 81 has the same conditions as Runs 87, 88, 90, and 91 except it is conducted with unsieved JSC-1 A. 
By fitting Equation (12) to the data, we can get a value for the factor F for unsieved JSC-1 A at 1000 °C 
which can be related to the effective particle size of the unsieved simulant. The value of F is found to be 

F= (1/90,000) 5 _1 for Run 81 


Effects of Temperature and Velocity 

Runs 8 1 , 80, and 82 have been conducted with unsieved simulant but with different temperatures, 
that is, 1000, 900, and 800 °C, respectively, with a flow rate of 0.5 cm/s (evaluated at standard conditions). 
Runs 74, 76, and 77 are duplicates of the previous runs except that the flow velocity is increased to 1 cm/s 
(evaluated at standard conditions). Figures 5 and 6 plot the conversion as a function of time for two of the 
temperatures (900 and 800 °C, respectively) and for both the 0.5 cm/s and 1.0 cm/s cases. Also plotted in 
the figures are the corresponding fits from the model. The following values of F are obtained: 

F= (1/90,000) s" 1 for 1000 °C 
(1/150,000) s' 1 for 900 °C 
(1/430,000) for 800 °C. 


+ 

Run 76: 

: Icm/s 

V 

Run 80: 

0.5 cm/s 

o 

Model: 

Diffusion 


0.7 - 

c 0-6- 


O 

V) 



Time.hr 


Figure 5. — Measured and model predicted conversion of unsieved 
JSC-1 A particles at 900 °C for two different hydrogen flow rates. 


NASA/TM— 2009-21 5618 


9 


1 


0.9 - 


0 Run 77:1 cm/s 
V Run 82: 0.5 cm/s 
O Model: Diffusion 


0.8 - 


0.7 - 


c 06 - 
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Figure 6. — Measured and model predicted conversion of unsieved 
JSC-1 A particles at 800 °C for two different hydrogen flow rates. 


It is clear from the data that as the operating temperature is reduced, the conversion in a given amount of 
time is also reduced. Initially, the higher flow rate results in a higher conversion rate. However, after 1 to 
2 hr, the higher flow rate and lower flow rates yield essentially the same conversion in agreement with the 
model. 


An Arrhenius fit for F, that is, a linear relationship between log(F) and \!T where T is the operating 
temperature in degrees Kelvin, may be derived from the PNNL data. One of the reasons for employing an 
Arrhenius fit is that the factor F contains the equilibrium constant k which has in the past (Ref. 10) been 
fitted by an Arrhenius expression. The equation for F is found to be 


and it is plotted in Figure 7. An approximate error band (±7 percent) in log(F) based upon the spread in 
the data is also shown in Figure 7. This is equivalent to approximately ±20 percent error in F. 


The LM data have been obtained in a much larger scale setup than the PNNL data, as noted earlier. 
However, all of the factor parameters that control the factor F are the same for both setups. The only 
exception is that some of the LM tests have been conducted at temperatures lower than 800 °C which was 
the lowest temperature for the PNNL tests. The equation for the factor F (Eq. (13)) is strictly valid in the 
range of 800 to 1000 °C which is the temperature range of the PNNL experiments. The LM experiments 
which are in this temperature range are categorized as in-range and those at lower temperatures are 
categorized as out-of-range. Both categories are considered in the following discussion. 


Fit for F 


log (F)=- 



(13) 


T 


LM Data and Model Validation 
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Figure 7. — Arrhenius fit for the factor F for unsieved JSC-1 A in the 
temperature range of 800 to 1000 °C (i.e., 1073 to 1273 K). 



Figure 8. — Measured and predicted conversion versus time for 
unsieved JSC-1 A at 970 °C. 


In-Range Predictions 

In-range LM data are available at temperatures of 970, 950, 925, 900, and 800 °C. Figures 8 to 12 
compare the model predictions of conversion versus time with the measured data. As before, a conversion 
of unity represents 2.8 percent water yield. Predictions are made using values of F at the higher and lower 
error bars. With the exception of some discrepancy at later times for the 900 and 925 °C cases, the 
“in-range” LM data are in very good agreement with the model predictions. It may also be noted that the 
predictions at earlier times in the conversion process appear to agree better with the LM data than with 
the PNNL data. This may be due to better mixing of hydrogen with the simulant in the rotating reactor 
configuration as compared to the tube reactor which may lead to decreased effects of residence time 
(see Eqs. (10) to (12)). 
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Figure 9. — Measured and predicted conversion versus time for 
unsieved JSC-1 A at 950 °C. 


1 

0.9 

0.8 

0.7 

c O.B 
Q 

'« 

® 0.5 
<3 0.4 
0.3 
0.2 
0.1 


I = 

■■■Model: High 
■■■Model: Low 
▼ Measured 


— 


5>V* 







T ▼ 


▼ T 


0.5 


1.5 


2.5 


3.5 


Time, hr 


Figure 10. — Measured and predicted conversion versus time for 
unsieved JSC-1 A at 925 °C. 
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Figure 11. — Measured and predicted conversion versus time for 
unsieved JSC-1 A at 900 °C. 



Figure 12. — Measured and predicted conversion versus time for 
unsieved JSC-1 A at 800 °C. 
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Out-of-Range Predictions 

Out-of-range LM data are available for temperatures of 700 and 600 °C. For lack of any alternatives, 
the Arrhenius fit for F (Eq. (13)) is used. Model comparisons with the LM data are shown in Figures 13 
and 14. It is seen that the agreement remains good at 700 °C. However, at 600 °C, the model with the 
assumed fit for F underpredicts the results. It is possible that there may be temperature variations in the 
reaction rate and F which are not captured by the model or the Arrhenius fit at this temperature. Note that 
experiments show that 600 °C appears to be a lower temperature limit for measurable hydrogen reduction 
of JSC-1A (Refs. 1 and 16). 
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Figure 13. — Measured and predicted conversion versus time for unsieved 
JSC-1 A at 700 °C. 
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Figure 14. — Measured and predicted conversion versus time for unsieved JSC-1 A 
at 600 °C. 
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Conclusion 


A model has been developed for hydrogen reduction of JSC-1 A and validated against data from two 
separate investigations with vastly different experimental scales. It is shown that a single factor, F, related 
to the gas diffusion coefficient inside the particle, controls the conversion rate at high conversion fraction 
or for small gas residence time in the regolith bed. Decreasing the particle size and increasing the reaction 
temperature both result in more rapid conversion under these conditions. An Arrhenius fit for the factor F 
has been obtained and shown to be accurate in the temperature range of 700 to 1000 °C (i.e., 973 to 
1273 K). The developed model can serve as a starting point for understanding the hydrogen reduction 
process for various types of lunar regolith but experimental data will be required to obtain the factor F 
for these regolith types. This model is also useful for overall hydrogen reduction system studies and 
optimization. 
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